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A general orbit-averaged guiding-center Fokker-Planck operator suitable for the numerical analysis 
of transport processes in axisymmetric magnetized plasmas is presented. The orbit-averaged guiding- 
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I. INTRODUCTION 

Transport processes play a dominant role in the long-time behavior of strongly-magnetized plasmas. In the absence 
of wave-induced (or turbulent) transport, the long-time magnetic confinement of plasmas is based on the small 
dimensionless parameter es ^ p/ Lb ^ 1 (the ratio of the characteristic gyroradius p and the magnetic nonuniformity 
length scale Lb)- For such plasmas, the dimensionless parameter = Lb/\v {^u being the mean- free-path) can be 
used to describe different classes of coUisional transport processes, such as classical coUisional transport p| {Xi, <C 
Lb) and neoclassical "collisionless" transport [2] (A^, ^ Lb)- The quasilinear (wave- induced) transport processes 
associated with rf-induced heating and current drive, on the other hand, cause a slow time evolution of the background 
plasma distribution (as a result of one or more wave-particle resonances) with a time-scale ordering that is quadratic 
in a small parameter — |Buj|/|Bo| <C 1 associated with the rf-wave amplitude [3|,[1]- 

Because the long-time behavior of the plasma distribution function for each particle species depends on competing 
coUisional and quasilinear transport processes, an accurate treatment of both transport processes in realistic magnetic 
geometry is a crucial element in determining the equilibrium and behavior of fusion plasmas. For this purpose, the use 
of dynamical-reduction methods (e.g., guiding-center transformation 0,01) can yield reduced transport operators (in 
phase spaceVthat possess attractive numerical properties in addition to accurately representing coUisional and/or 
quasilinear *4"| transport processes of interest. 

The purpose of the present paper is to present a brief derivation of a general orbit-averaged guiding-center Fokker- 
Planck operator suitable for numerical studies of transport processes in general axisymmetric magnetic geometry. 
This reduced Fokker-Planck operator represents drag and diffusion processes in a three-dimensional space composed 
of guiding-center invariants. In this reduced formulation, the conjugate orbital angles have either been eliminated 
from the guiding-center Fokker-Planck operator by averaging or are absent by axisymmetry. 

The remainder of the paper is organized as follows. In Sec.|lll we discuss the guiding-center Hamiltonian dynamics 
of charged particles in unperturbed axisymmetric magnetic geometry. In Sec. IIIIl we introduce the guiding-center 
Fokker-Planck operator previously derived for general magnetic geometry and arbitrary guiding-center orbit topology 
[2]. In Sec. IIVI we first introduce the orbit-averaging operation for standard and non-standard guiding-center orbits 
in axisymmetric tokamak geometry. Next, we present the orbit-averaged guiding-center Fokker-Planck operator and 
discuss its properties. In Sec.|Vl we briefly discuss the derivation of the bounce-center Fokker-Planck operator obtained 
by performing the bounce-center phase-space transformation Q on the guiding-center Fokker-Planck operator and 
discuss its connection to the orbit-averaged guiding-center Fokker-Planck operator derived in the previous section. 
Lastly, we summarize our work in Sec. IVII and discuss its applications. 



II. GUIDING-CENTER DYNAMICS IN AXISYMMETRIC MAGNETIC GEOMETRY 

The existence of the small parameter ee ^ 1 in magnetically-confined plasmas forms the basis of the unperturbed 
guiding-center dynamical reduction [6|, in which the fast gyromotion time scale associated with the gyroangle Cg 
(with C,g = eg^ri), with the gyroaction Jg = fiB/V, acting as its canonically-conjugate (adiabatic) invariant, is 
asymptotically decoupled from the parallel and cross-field motions of a guiding-center particle. The unperturbed 
guiding-center dynamics is expressed in terms of the guiding-center position X and the parallel velocity U|| , where the 
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guiding-center velocity 



^ b / ^ 

X = W|| b + ee — — X (^/i VB + mi;| b • Vb ) = b + eb vb (1) 



is decomposed in terms of a parallel velocity v^^ = b • X along a field line and a slower cross-field drift velocity Vb due 
to weak magnetic-field nonuniformity. The guiding-center parallel acceleration 



f II = — 



m \ till 



EB — • VB (2) 



guarantees that the guiding-center kinetic energy 

is a constant of motion {£ = for a time-independent magnetic field and in the absence of an electric field) . In Eqs. 
0, the guiding-center magnetic-moment invariant /i and the coordinates (X, in the reduced four-dimensional 
guiding-center phase space are expressed as asymptotic expansions in powers of the small parameter eb , where first- 
order corrections explicitly take into account magnetic- field nonuniformity Q. 

We note that, while a nonrelativistic guiding-center formulation is considered here, its generalization to a relativistic 
formulation (appropriate for fast electrons) can easily be accommodated fld\. The kinetic energy ([3]) is thus replaced 
with 5 = (7 — l)mc^, the parallel velocity is replaced with the relativistic parallel momentum p|| — jmv\\, and 
the magnetic moment /x is replaced with the relativistic magnetic moment /i = \p±\'^ /2mB, where the relativistic 
factor is 7 = (1-1-2 fMB/mc^+pj /m^c^Y^'^. Additional details on relativistic guiding-center dynamics in axisymmetric 
magnetic geometry can be found in Ref. [Tlj . 

A. Axisymmetric Magnetic geometry 

The general a xisy mmetric magnetic field considered in the present paper is expressed in terms of three equivalent 
representations [12| : 

X vip + q{tp) vip X ve 

B = { B^\/(j) + BgS/9 + B^VV' (4) 
B't' d^/d(t) + B<^ dX/de 

where the two-covariant, covariant, and contravariant representations (from top to bottom, respectively) are expressed 
in terms of the (poloidal) magnetic flux ip, which satisfies the condition B • V?A — (i.e., magnetic field lines lie entirely 
on a constant-?/' surface), and the poloidal and toroidal angles 6 and 4>. The toroidal and poloidal components of the 
magnetic field Q are B^^, = B^/R = B'f' R, where R = |5X/90| = |V0|-\ and Bpoi = \Vi^\/R = B^ \dX/de\. The 
safety factor q{^) appearing in the two-covariant representation is defined as [isj 

_ B-Wcf, B^ 
"^"^^ = bTw = B^- 

We note that, because of the magnetic-fiux condition B • Vf/; = 0, the covariant component B^ = 
— Bg {Wd in Eq. ^ vanishes only if the coordinates ip and 6 are orthogonal. The spatial Jacobian 

V associated with the coordinates {ip,6,ip) is 

V = {Wtl^xVe-V4>)-^ = (B-V6l)-i = (B^)"\ (6) 

where B^ is assumed to be positive. Lastly, the infinitesimal length element along a magnetic field line is 

ds ^ ^ d0. (7) 

Note that, according to the standard axisymmetric tokamak ordering (with B ~ i?tor), the ratio B/B^ can also be 
expressed as B/B^ — B / [B tor /{qR)] — qR, so that we recover the standard approximation ds ~ qRdO. 
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B. Guiding-center Motion in Axisymmetric Magnetic geometry 



The guiding-center motion in arbitrary magnetic geometry described by Eqs. HI)-© possesses two constants of the 
motion: the total guiding-center energy ([3]) and the guiding-center magnetic moment fj,. Guiding-center motion in 
general axisymmetric magnetic geometry ^ is also characterized, according to Noether's theorem, by a third constant 
of the motion: the toroidal canonical guiding-center momentum 

Here, the vector potential A — + ■(/'tor(V') (with q = di/jtov/dip defined in terms of the toroidal magnetic 
flux tptoi) was obtained from the two-covariant representation in Eq. ([¥]), where the coordinates {ip, 9, (p) now describe 
the guiding-center position X, py = v\\/Q. denotes the parallel gyroradius, and -B<^ = B • dl^/d4> denotes the covariant 
component of the axisymmetric magnetic field. 

The projection of the two-dimensional drift surface 

^ - p,l i?^ = (9) 

onto the poloidal plane (X, Z) generates a closed curve ip — tp^O) parameterized by the poloidal angle 9 and labeled 
by the guiding-center invariants (■!/), £,/i): 

ipie,a;lp,£,fi) = lp + S^i9,a;lp,£,fi), (10) 

where the "bounce-radius" Sip = tJj — ^ represents the departure of the drift surface (jlOp . labeled by ^, from a 
magnetic surface "0 and a — ±1 denotes the sign of . A guiding-center orbit O is obtained either by integrating the 
guiding-center equations of motion 



?/; = X • Vip = W|| b • + ee vb • ^ip = cb "Pb 

9 = ±-V9 = v\\h-V9 + eBVB-V6l = v\\ByB + eB ^b 



(11) 



for a given set of guiding-center invariants {ipjE, ^) or generating the orbit directly by the constant-of- motion method 
[l^ . For each generic set of guiding-center invariants there corresponds a unique guiding-center orbit O, 

which is either a trapped-particle orbit (if the bounce-radius Sip vanishes along the orbit), or a passing-particle orbit 
(if the bounce-radius 5ip does not vanish). Non-generic orbits [T"5i] include the stagnation orbits (where ^ = or 
v\\B^ / B — — eB 9b 7^ 0) and the pinch orbits (barely-trapped orbits in the zero-banana- width limit). 

On each flux-surface ip, the magnetic field amplitude is assumed to vary monotonically between a minimum value 
Bo{ip) and a maximum value Bi{ip), located at the poloidal locations 9(j{ip) and 6*1(7/)), respectively. The turning 
points (V'b, ) of ^ trapped-particle orbits {Tp,£,fi) are located on the drift-surface label iph = ip, where the magnetic 
field reaches its maximum value along the orbit and S'ip{9^ ^ili'ipi.jS , fx) = 0. Consequently, the X point separating 
trapped-particle and passing-particle orbits parameterized by the drift-surface label ip is located at the position 
[V^,^i(?)]. 

While the magnetic moment fi is an important invariant for guiding-center dynamics, the pitch-angle coordinate 



C(V^;£,m) = v\\/v = cj sjl - tiB{iP)/£. (12) 

is better suited in describing the transition between trapped-particle and passing-particle orbits. In order to convert 
the pitch-angle coordinate into a suitable guiding-center invariant, we replace the guiding-center magnetic-moment 
invariant /i with [16] 



= •\/l-/iSo(^)/5. (13) 

The physical interpretation of the pitch-angle invariant ^0 can be given in terms of its connection with the bounce- 
action invariant Jb. With this definition, the trapped-passing separatrix given by the relation, 

1 - [l-il)BM/BM^^ 



which is an even function of ^0 a-nd does not depend upon the energy £ (important for numerical applications). 
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The guiding-center Jacobian in coordinates (X,p, ^, Q) is J7gc = p^, while in terms of the coordinates (X,p, ^o, Cg): 
it is Jgco = \d^/d^Q\ —p^ *Co/|'?|, where we used the definition ^(^,0) = 0)/-Bo(V') with the relation 

lei - .^1 - ^ (14) 

Next, the spatial Jacobian in coordinates X = ("0, 0, 0) is V = [B^ {ij])]^^ , so that the total guiding-center Jacobian in 
coordinates (%l^,9,(f>; p, <^o,Cg) is 

In the next section, we express a general guiding-center Fokker-Planck operator in terms of the guiding-center invari- 
ants ("07^1 Co) in a- form that will be suitable for orbit averaging in Sec. IIVI 



III. GENERAL GUIDING-CENTER FOKKER-PLANCK EQUATION 



The Fokker-Planck equation 0,11 [T3| forms a paradigm for the investigation of classical, neoclassical, and quasilincar 
transport processes in plasmas. When written in terms of the test-particle phase-space coordinates (x, p), the general 
Fokker-Planck operator is expressed as a local partial-differential operator in momentum space [2] : 



C[/](x,p) = 



-|-(K/(x,p) 



D- 



g/(x,p) 
9p 



(16) 



where the Fokker-Planck momentum- friction vector K and the momentum-diffusion tensor D are functions of (x, p) . 
When representing coUisional transport processes, for example, these Fokker-Planck coefficients K = K[/'] and 
D = Y^' D[/'] are expressed as integral operators acting on the field-particle distribution /' (where the field-particle 
species may coincide with the test-particle species). Hence, the Fokker-Planck operator (|16p may either be a linear, 
bilinear, or nonlinear operator acting /, depending on the type of transport problem one wishes to investigate. 

While collisions and wave-particle interactions take place locally in particle phase space z — (x, p) , the transfor- 
mation to reduced phase-space coordinates z (e.g., guiding-center coordinates) will generically introduce transport 
coefficients in the full reduced phase space. This implies that the Fokker-Planck coefficients (if* , D'-' ) in three- 
dimensional particle-momentum space are replaced with six-dimensional reduced phase-space Fokker-Planck coeffi- 

defined as 



cients (a , _D ) 



K (z) 



A-(z) 



az"(z) 



and 



=z(z) 

and the Fokker-Planck operator p6)) transforms to 



9p* 



dpi 



z(z) 



C[/](z) = - = 



1 d 



(17) 



(18) 



where /(z) = /(z) is the particle distribution expressed in terms of the reduced phase-space coordinates and J is the 
Jacobian for the transformation z — s- z. The most important aspect of the transformation to Eq. (|18p involves the 
choice of the new phase-space coordinates z. Note that the definitions (fT7| for the reduced phase-space Fokker-Planck 
coefficients require that the transformation z z(z) and its inverse z ^ z(z) must be known up to the desired order 
in magnetic-field nonuniformity. 

One possible choice is to adopt a canonical action-angle formulation [l^, where z = (J,0) includes the three- 
dimensional action invariants J for magnetically-confined particles and their canonically-conjugate angles (which 
are ignorable coordinates by construction df/dO = 0). While the three-dimensional action-space Fokker-Planck 
operator derived by Bernstein and Molvig 18] formally describes classical and neoclassical transport processes in 
axisymmetric magnetic geometry, it is not suitable for numerical implementation since some of the action coordinates 
(e.g., the bounce action Jb) are not local coordinates. Hence, the action-space Fokker-Planck operator does not 
describe local transport processes, which makes the transport analysis of its results difficult to interpret. 

Another choice is to adopt local noncanonical guiding-center coordinates 0, [lB| leading to the construction of a 
reduced guiding-center Fokker-Planck operator. In a weakly-nonuniform plasma with a strong axisymmetric mag- 
netic field, for example, the reduced guiding-center Fokker-Planck equation describes transport processes in a four- 
dimensional space: the poloidal-flux and poloidal-angle coordinates (?A,^) in physical space, and the energy and mag- 
netic moment coordinates (£, fi) in velocity space, which are invariants for the guiding-center motion. Here, the reduc- 
tion from six to four dimensions is associated with the fact that the guiding-center plasma distribution is independent 
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of the gyroangle (by definition) and the toroidal angle (f> (by axisymmetry) , so that transport along these ignorable 
angles is irrelevant. While Zaitsev et al. [T6| considered the lowest-order definitions (in magnetic- field nonunifor- 
mity) for their choice of guiding-center coordinates, Brizard considered first-order corrections (ee = p/Lb <C f) 
as well. In both works, classical collisional transport appears in Eq. ([TT]) as a result of the lowest-order relation 
i9X/9p = — b X I/mfi for the guiding-center position X(x, p) = x — Po(^iP)j where Pn = b X p/mfi represents 
the lowest-order gyroradius. The higher-order guiding-center corrections kept by Brizard [3] are consistent with the 
low-collisionality regime, where the mean-free-path > Lb is longer than the magnetic-nonuniformity length scale 
Lb. 



A. Guiding-center Fokker-Planck Operator 



The derivation of a reduced Fokker-Planck operator that is suitable for numerical implementation must begin with 
finding local invariant coordinates that also allow relative computational simplicity for realistic magnetic geometries. 
In the present work, the exact invariants associated with the time-independent axisymmetric magnetic geometry are 
the kinetic momentum p = and the drift-surface label ip = tph (for trapped particles) or tpt (for passing 

particles). As our third invariant coordinate, we use the pitch-angle invariant defined in Eq. (|13p . which allows an 
explicit representation of the trapping and detrapping transport processes. 

Once a set of invariant guiding-center coordinates /° — {ip,p,^o) is chosen (not necessarily action coordinates), the 
transformation of the Fokker-Planck operator can be greatly simplified by writing it in Poisson-bracket form as 



C[f] - ^ - {x\ (k^ f - D^' {x\ /})}, (19) 

where the noncanonical Poisson bracket { , } is used to replace momentum partial derivatives dg/dpi = {x^,g}. 
The significant computational advantage of this Poisson-bracket formulation (fT51) is based on the fact that Poisson 
brackets transform naturally under the type of phase-space transformations considered in the present work (i.e., those 
generated by Lie-transform methods). 

The general guiding-center Fokker-Planck operator was derived by Brizard Q as a result of the guiding-center 
dynamical reduction of the Fokker-Planck operator ((TF 



where Tgc and T^^^ denote the pull-back and push-forward (Lie-transform) operators associated with the guiding- 
center transformation for a nonuniform magnetic field, (• • •)g denotes an average with respect to the guiding-center 
gyroangle Cg, and { , }gc denotes the guiding-center Poisson bracket. In Eq. ([20|) . the generalized gyroradius vector 

contains first-order corrections associated with magnetic-field nonuniformity. 

When we use the phase-space-divergence property of Poisson bracket 



G 

' gc 



where F and G are two arbitrary guiding-center phase-space functions, we can write the divergence form of the 
guiding-center Fokker-Planck operator ([20|) : 



C,c[F] E ^ ^ 



J ( /Cg", F ~ Pgf ^ 



(21) 



JdZ 

where J' is the total Jacobian (|15p . The guiding-center Fokker-Planck friction components 

^gc ^ (Tg/K.A")^ (22) 
and the guiding-center Fokker-Planck diffusion components 

Pgf EE /(A")^.Tg-iD.A^\ , (23) 



are expressed in terms of the guiding-center push- forward expressions Tgj.^K and Tgj.^D of the particle momentum- 
space Fokker-Planck friction vector K and diffusion tensor D, and the vector- valued projection coefficients 

A" = {X + p„ ^ A^.VZ" A^ ^ + A^ — . (24) 
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Here, the guiding-center vector-valued projection coefficients in (X,£,/i) space Q 

A' = pjm \ (25) 

= (n/B) dpJdQ J 

are constructed from the guiding-center Poisson-bracket expressions involving the guiding-center push-forward of the 
particle position T^^^x = X + p^, where the generalized gyroradius = T^^p and the generalized guiding-center 
momentum 

contain first-order corrections associated with the guiding-center transformation. Note that, while the projection 
coefficients A'^ and A'' in Eq. (|25p retain first-order corrections in magnetic-field nonuniformity, we have omitted 
first-order corrections in A''^ since they yield second-order contributions in the final expressions for Eq. (|2ip . 



B. Guiding-center Projection Vectors 

Using the definitions for the invariant-space coordinates = ('0,p, Co)i the guiding-center vector- valued projection 
coefficients A° = (A'^, A^, A*°) are expressed as 



A"^ = A^-Vt/; = — X Vt/;, 
mil 



AP 
A& 



Pe 



A^.V^o + A^?| + A^^ = 



5^0 



2^0 



Pe ^ dPc 



din Bp 
(hp 



(27) 
(28) 

(29) 



We now discuss the physical nature of the projection vectors ([27|)-(f29l 



1. Guiding- center Radial Projection 
Using the magnetic representation ([5]), the guiding-center radial projection vector (j27p becomes 

where we introduced the angle-like coordinate x = — q{ip) 0, such that the two-covariant representation of the 
axisymmetric magnetic field ^ becomes B = Vx X Wip. Here, = (Vi/' X V^) X V^/'/l Vi/'P denotes the projection 
of the gradient of 6 that is perpendicular to Vip, while = V0 (since V(p-Vip = 0). Hence, the guiding-center 
radial projection vector ([5U)) projects momentum-space transport processes onto a magnetic-flux surface in a direction 
that is locally perpendicular to b. 

The guiding-center radial projection vector A''', whose magnitude is 

_ \S/lp\ cBpoi 

therefore generates the local projection of the guiding-center push-forwards of the momentum-space Fokker-Planck 
friction 

4 = (Tg-c^K)^ • X (31) 
and the momentum-space Fokker-Planck diffusion 

V^J = {A'^ . T-iQ)^ ~ X VV^, (32) 



\ SC /g j^Q 
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where a ^ ip va Eq. p2|) . These projections only take into account transport processes occuring within a magnetic- 
flux surface. The connection between radial transport and toroidal and/or poloidal rotation is therefore naturally 
contained within this projection. 



2. Guiding- center Energy Projection 
The guiding-center energy projection vector (|28p is expressed as 

AP ^ -P^ = p + eAf + (34) 

where denotes the first-order correction of A^. It can be explicitly expressed in terms of the gyroangle-independent 
first-order gyroradius Pg^. = (b/il) X vb and = + e as 

where p = p/{mD,) is the gyroradius magnitude for a deeply-trapped particle (i.e., for = 0) and 

i (l-e) VlnS + b-Vb . 

It can be shown explicitly, however, that Aj* does not have any component directed along the zeroth-order unit vector 
p (because the guiding-center kinetic energy is identical to the particle kinetic energy). Hence, the magnitude of the 
guiding-center energy projection vector (p8)) is | A^j = 1 up to second order in magnetic-field nonuniformity. 



p X b = phx 



3. Guiding- center Pitch-angle Projection 



The guiding-center pitch-angle projection vector (|29|) is expressed as 



%^-fbxVlni.o 
Po (1 - Co) 2 



_ - Co C C , ^ A Co 



(36) 



where A^" denotes the first-order correction of A^". To lowest order in magnetic-field nonuniformity, we find |A^°| — 

yT^C/(pV^Co). 



C. Guiding-center Fokker-Planck Equation 



Now that the guiding-center Fokker-Planck components are expressed in terms of the invariant-space coordinates 
= {'4'^P^ Co) I we turn our attention to the guiding-center Fokker-Planck kinetic equation. When the original guiding- 
center coordinates (V', 9, 0; f , /i, Cg) are used in describing guiding-center dynamics in axisymmetric magnetic geometry 
(where we replace the guiding-center parallel velocity ?;|| with the guiding-center energy f and the gyroangle C^^ is an 
ignorable coordinate [3), the guiding-center Vlasov evolution operator is ordered as Q 

-dt ='^d^^ P -B + 89 + ^^^^ 9^' (^^^ 

where we used d/dcj) = (by definition of axisymmetry) , explicit time dependence is over long time scales (i.e., 
d/dt = erd/dr is ordered small), and the guiding-center dynamics on the poloidal plane is represented by Eq. pT|) . 

For each generic set of guiding-center invariants (ip^S,^), the guiding-center Fokker-Planck kinetic equation is 
expressed as 



(38) 
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where the magnetic-flux drift motion 



is expressed in terms of the bounce radius ^TU\i when it is projected onto the guiding-center orbit O associated with the 
guiding-center invariants {tp,£,fj,). While the topology of these orbits can be rather comphcated [l4,, J5, 20, J.l,_22], our 
discussion will remain as general as possible concerning the nature of the axisymmetric magnetic geometry. Whether 
a guiding-center orbit O corresponds to a trapped-particle orbit or a passing-particle orbit, however, the motion is 
periodic in the poloidal plane {ip,0) and a generic orbital period in the poloidal plane 

f de 

TO = f ^ (40) 
Jo p 

is defined as a closed-loop integral along the drift-orbit curve O parameterized by the poloidal angle 9 (at constant 
invariants ip^ and fi). Note that the definition (|40)) of the orbital period can also be given in terms of dip/ip 
[31 which, when evaluated along the guiding-center orbit O, yields {d^p/tp)o = d9 {dS'ip/d9)/ip — d9/9 upon using 
Eq. dSSl). 



IV. ORBIT-AVERAGED GUIDING-CENTER FOKKER-PLANCK EQUATION IN 

LOW-COLLISIONALITY REGIME 

A. Orbit- Averaging Operator in Axisymmetric Magnetic Geometry and Low Collisionality Regime 

In the neoclassical transport regime e^, — vto ^ 1, the guiding-center Fokker-Planck kinetic equation (j38p reduces 
to leading order to 9 dF/d9 = 0, which states that F = F{'tp,£ , fi;T) is independent of the poloidal angle 9. The 
physical picture is that guiding centers orbits undergo many poloidal cycles before being perturbed by collisions. 

We now introduce the or&ii-averaging operation 

{■■■)o ^ - i (--of , (41) 
To Jo 9 



where the orbital period to is defined in Eq. ()40p . By orbit-averaging the guiding-center Fokker-Planck kinetic equation 
(I38p in the low-coUisionality approximation, we finally obtain the orbit-averaged guiding-center Fokker-Planck kinetic 
equation 

^r-Q^ = (Cgc[i^])o, (42) 
which describes the collisional time evolution of the orbit-averaged guiding-center distribution 

{F)o - F(^,£,M;r), (43) 



and the guiding-center Fokker-Planck operator is given by Eq. ([2T|) . 

Lastly, we note that the relation between the magnetic flux ■0 and its orbit-averaged value {ip)o (the "drift-center" 
position) is expressed as 

^ = (^)o + (S^ - {H)o) , (44) 

where the drift-surface label is '0 = {ip)o — (<5V')o- For trapped and passing particles, we therefore obtain the following 
relations for the deviations from a magnetic-flux surface: 

( 5i! \ 

(45) 

where we used the fact that ((50)b = for trapped particles, while {54')t 7^ for passing particles (since ^0 does not 
change sign along a passing-particle orbit). These relations show that, while the deviation for a trapped particle 
may be large, the deviation (5^j — {5ip)t) for a passing particle is in general small. 
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B. Approximate Orbit- Averaging Operator 



In this Section, an explicit expression for the orbit-averaged guiding-center Fokker-Planck operator that appears on 
the right side of Eq. (|42p is derived. In general, the orbit-averaging operation (j4ip must be computed along numerically 
calculated orbits. However, useful analytical expressions of the orbit-averaged guiding-center Fokker-Planck operator 
can be obtained when the orbit deviation (5'0/|V'(/'| is small as compared to the local minor radius r. The difference 
between the exact incremental time element dO/O and the approximate time element ds/wy is of order e| and can be 
neglected when compared to the corrections of order ee in the Fokkcr-planck operator. The orbital period (j40[) thus 
becomes 

Jo n « Jo B%ij,e) 

The orbit topology is explicitly taken into account in Eq. (|46p through the dependence of the integrand B/{(_B^) on 
the bounce-radius 5ip{0,'j\%lj,£, ^) = V-" — defined in Eq. pop . where ^ = (for trapped particles) or tpt (for passing 
particles). Note that the orbital period (l46l) is exact for stagnation orbits (where ^ = yields = — I'll / B). 
In the definition of the orbital period (|46p, it is convenient to extract the magnetic- geometric factor 

A(^,p,Co) ^ r de^^ (47) 

which defines the length of an orbit on the drift-surface labeled hy ip. In the standard axisymmetric tokamak ordering, 
we find A ~ 2'k qlip) Rq, where i?o denotes the major radius of the magnetic axis. Next, we define the normalized 
orbital period 

X{^,P,Co) ^T<f del^^(^^)ro, (48) 



A Jo ^B<> \ A 

which satisfies the condition X{ipt,P, I'Col = 1) = 1 for completely-passing particles (i.e., |^| = 1 and /i = so that 
= V is a constant of motion). Note that the orbital period pS)) [or normalized orbital period (pS)) ] becomes infinite 
on the boundary that separates trapped-particle orbits and passing-particle orbits (i.e., pinch orbits). Hence, the low- 
collisionality approximation {dF/dO — 0) does not technically hold very close to the trapped-passing boundary, where 
the bounce and transit periods become much larger than the characteristic collisional time scale . However, the 
fraction of particles with A ^ 1 is very small since the corresponding singularity is integrable. Indeed, by introducing 
^ob such that 1 - * (1 - ^^^J ^ 0, it can be shown that 

A(^,P,^o)d?o = ^ j ^ ^dO J^^^ = J j ^ ^dO^ 

which is a geometrical factor of order unity. 

Using the Jacobian (fT5|) . the normalized bounce period ((48|) is now expressed as 

A(^'P'^o) = ^ (f Jdd = 2'KJo (49) 

so that the orbital period (l46l) becomes 



2tiJoBq 

TO = 57 — . (50) 



The orbit-averaging operation (j4ip becomes 



{■■■)oi^,P,^o) ^ ^ i (•••)- = - / (•••)4S - ^ /(•••)^?, (51) 



TO Jo TO Jo v^B'^ Jo Jo 27r 

For trapped-particle orbits, the orbit-average (jSip yields the explicit formula for bounce-averaging operation 

'^b ^ Jo7 «KI ^ Jo ^ ^ Je- 27r 
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Here, the summation is over the sign cr of ^ and, therefore, the symbol ^ takes the average (for trapped particles) 
between values of the integrand for ^ > and values of the integrand for ^ < 0. For passing-particle orbits, on the 
other hand, the orbit-average (|5ip yields the explicit formula for transit-averaging operation 



1 f^"' , . ^ de 



C. Orbit-averaged Guiding-center Fokker-Planck Operator 



In the expression for the guiding-center Fokker-Planck operator (|2ip , the guiding-center function F is independent 
of the angle coordinates {9,(j),Cg)i hi the limit of weak coUisionality and under the assumptions of axisymmctry and 
gyroangle invariance. Hence, only derivatives of the guiding-center distribution F with respect to the guiding-center 
invariants /° will remain in Eq. pTjl . The coefhcients (A^gci^^g^) and the Jacobian J' in the guiding-center Fokker- 
Planck operator (|2ip . on the other hand, still depend on the poloidal angle 9. We therefore need to orbit-average the 
guiding-center Fokker-Planck operator (PT|) according to the orbit-averaging procedure (|5ip . 

The orbit-averaged guiding-center Fokker-Planck operator is expressed as 



dF 



Jo {^lc)o F - (^gc)o 7m 



(54) 



where the averaged Jacobian J^q now becomes the Jacobian for the orbit-averaged guiding-center Fokker-Planck oper- 
ator (|54p . Orbit-averaged guiding-center Fokker-Planck operator that treat friction and diffusion in three-dimensional 
space have been derived previously p^ . While these previous operators retained only the lowest-order terms in 
magnetic-field nonuniformity, the orbit-averaged guiding-center Fokker-Planck operator (|54p retains first-order cor- 
rections as well, which allows a realistic magnetic geometry to be considered. 

The orbit-averaged guiding center Fokker-Planck (FP) operator is widely used in the physics of magnetized plasmas 
(earth magnetosphere thermonuclear fusion [l3,[l3[, .. ). For practical applications (hke studying fast particle 
dynamics generated by rf waves or a constant electric field [24|). analytical expressions often correspond to an over- 
simplified description of the problem of interest, and therefore a more realistic approach requires full numerical 
calculations. 

For computational purposes, the orbit-averaged guiding-center Fokker-Planck operator (jM]) is recast in a flux- 
conservative form 



{CAF])o ^ 



J_d_ 

Jo dTp 



{jo\v^\os;^ 



]_d_ 

2 dp 



{p' Si) 



«0 



(55) 



such that the usual two-grid discretization technique may be applied for the finite-difference method [25| . Particle 
conservation is therefore naturally satisfied numerically up to second order in the truncation error. In Eq. (j55p . 

describes particle flux across magnetic flux surfaces, while and S^" account for momentum and pitch-angle 
dynamics respectively. From Eq. ([5^ . we have 



Si 

\ St J 



( Kt\ 

Kl 

V ) 



/ j^iH, ^^p jy^i \ 



F - 



JJPi' JJPP J^pi 



d/dp 



(56) 



\-p-'V^^d/dCo J 



where the friction coefficients coeflacients and symmetric diffusion coefficients D^^ = D^°- may be expressed in 
terms of guiding-center friction and diffusion components ([22|) and ((23|) according to the general relations 



Kl 



Kt = - 



P 



, ,-Vln*-/C^ 



(57) 
(58) 
(59) 
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and 



Dl 



1 



1 



DfP = 



P 



2^0 \* 



VV'-X'^^-Vln* 



1 f f 



2^0 
e2 



(60) 

(61) 

(62) 
(63) 
(64) 

(65) 



The numerical implementation of the orbit-averaged guiding-center Fokker-Planck operator (|55|) is described in a 
forthcoming paper [25(] using a novel 3-D scheme with fully-implicit time evolution. It is incorporated in the code 
LUKE [2||, which was initially developed for fast electron physics, and will therefore extend its range of applicability 
to multispecies physics and transport processes. 



V. BOUNCE-CENTER FOKKER-PLANCK OPERATOR 



The next step in the derivation of a reduced Fokker-Planck operator is to proceed with the construction of a 
bounce-center Fokker-Planck operator 



CbclF] 



TbcF 



(66) 



where T^^ and Tbc are the push-forward and pull-back operators associated with the bounce-center phase-space 
transformation 6] (and references therein), F = T^^F denotes the bounce-center distribution (which is independent 
of the bounce angle Cb), and ( )b denotes averaging with respect to Cb- In the bounce-center Fokker-Planck operator 
([66|) . finite-orbit effects will explicitly be taken into account. 

To see how these finite-orbit effects might arise in the bounce-center Fokker-Planck operator ((66|) , we first express 
the bounce-center push- forward and pull-back operators as T^^ = exp(±eGi • d), where Gi • d = G° 9a is defined 
in terms of the components G° of the first-order generating vector field for the bounce-center transformation and 
e = ed — tij/T(\ ^ I denotes its ordering parameter. Next, we expand the transformed operator 



TbcF 



eGid 



gc 



F 



(67) 



in powers of e, where the bounce-center Jacobian is defined in terms of the guiding-center Jacobian J ci& J 
J — t da{Gi J) + • • •. By keeping terms up to second order in e, and rearranging terms, we obtain 



(j + e a,(G?^) + • • •) Cgc + e F 



F 



G\ dbF 



(68) 



where the first term is simply the original guiding-center Fokker-Planck operator multiplied by the guiding-center 
Jacobian J = J ~\- e da{GiJ) (i.e., the first term is an exact divergence), while the second term is also an exact 
divergence in bounce-center phase space. The total bounce-center Fokker-Planck operator is therefore guaranteed to 
retain its phase-space divergence form. 

In future work, it will be our purpose to show that finite-orbit effects will appear in the bounce-center Fokker-Planck 
operator (|66p in the form of the second-order term 



Gt'-l 

dip 



J / b 



where ?/; denotes the bounce-center position and G'f denotes the corresponding component of the bounce-angle- 
dependent bounce-radius. In the present work, {'4')o = V' plays the role of the bounce-center position ip for a 
trapped-particle orbit {5^p plays the role of G^)^ while the nonlocal bounce-action Jb (canonically conjugate to the 
bounce angle Cb) is represented by the minimum-B pitch-angle coordinate ^q. 
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VI. SUMMARY 

The present work presents the derivation of a general guiding-center Fokkcr-Planck equation in axisymmetric 
magnetic geometry psp that aUows first-order corrections in magnetic-field nonuniformity to be retained. The Fokker- 
Planck operator (|2ip is obtained from a guiding-center transformation using the set of invariants (Tp,p,£^o). Next, the 
orbit-averaging procedure is introduced and an explicit Fokker-Planck operator (|54p is derived in the low-coUisionality 
limit. This operator is also expressed in a conservative form (|55|) that is best suited for numerical application. 

In the continuation of this paper (based on work briefly outlined in Sec. |V| , future theoretical work will consider 
finite-orbit effects in a different way by deriving a general bounce-center Fokker-Planck operator using Lie-transform 
methods. 
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